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Abstract — In recent years, there has been an increasing de- 
mand for ubiquitous streaming lilie applications in data networlis. 
In this paper, we concentrate on NUM-based rate allocation 
for streaming applications with the so-called S-curve utility 
functions. Due to non-concavity of such utility functions, the 
underlying NUM problem would be non-convex for which dual 
methods might become quite useless. To tackle the non-convex 
problem, using elementary techniques we make the utility of 
the network concave, however this results in reverse-convex con- 
straints which make the problem non-convex. To deal with such a 
transformed NUM, we leverage Sequential Convex Programming 
(SCP) approach to approximate the non-convex problem by a 
series of convex ones. Based on this approach, we propose a 
distributed rate allocation algorithm and demonstrate that under 
mild conditions, it converges to a locally optimal solution of the 
original NUM. Numerical results validate the effectiveness, in 
terms of tractable convergence of the proposed rate allocation 
algorithm. 

I. Introduction 

With recent advances in networking technologies and video 
compression, there is an increasing demand for ubiquitous 
multimedia appUcations Uke live streaming, video gaming, 
video conferencing, and voice over IP. Multimedia applications 
are characterized by a multitude of QoS requirements includ- 
ing stringent bandwidth, delay, and delay jitter guarantees. The 
ever increasing demand for streaming traffic has attracted a 
lot of research interests to develop efficient mechanisms for 
resource allocation between competing multimedia sessions in 
a wide variety of networking scenarios fTl-ISI. 

In the course of the last decade, rate allocation has been 
widely addressed as the (usually distributed) solution to Net- 
work Utility Maximization (NUM), which has emerged as 
an analytical framework to understand and design existing 
network protocols (|4l-||5l. The goal of NUM is to maximize 
the aggregate utility of the users subject to operational and 
practical constraints of the network. In the basic form of 
NUM proposed in jj), the feasibility of rate allocation was 
accommodated by congestion in links. So far, a plethora of 
studies have concentrated on NUM-based rate allocation for 
services with elastic traffic such as traditional file transfer Due 
to strict concavity and differentiability of the utility function 
for elastic traffic, such NUMs are smooth and strictly convex 
and thus far have been efficiently solved using dual or primal- 
dual methods (see e.g. ^ and references therein.) 



In contrast, applications that carry inelastic traffic like 
audio/video streaming, can only tolerate a limited amount 
of packet delay or fluctuation in rate. Hence, they are in 
possession of non-concave and often non-differentiable utility 
functions |[T|, |j6). This results in a non-convex and usually 
non-smooth NUM for which dual/primal-dual methods might 
prove quite useless. 

There have been several works that have addressed non- 
convex NUM problems for resource allocation supporting 
inelastic services |6 |-| 12 1. Lee et al. |7 | outlined the possibility 
of divergence of dual methods for non-concave utilities and 
proposed a distributed "self -regulating" heuristic for rate con- 
trol of non-concave utilities, where some of the sources turn 
themselves off according to their local information. Hande et 
al. fS"! proposed necessary and sufficient conditions for canon- 
ical distributed algorithm to converge to global optimum in the 
presence of non-concave utilities. A centralized algorithm for 
non-convex NUM has been proposed in [9| in which sum-of- 
squares technique was applied to a polynomial approximation 
of the non-concave utility function. However, this centralized 
approach suffers from high order of complexity. In fS], the 
authors exerted a redefined variant of the non-concave utility 
function in a distributed flow control algorithm so that the 
network can achieve a utility-proportional fair rate allocation. 
Authors of 1 10 1 merged the utility-proportional theory with a 
stochastic optimization framework to propose a rate control 
algorithm for the mixture of elastic and inelastic traffic in 
wireless sensor networks. In |11|, the authors introduced a 
smooth utility function as an approximation to the ideal stair- 
case utility function for SVC-encoded streams and leveraged 
the utility-proportional approach to redefine the NUM which 
is solvable using dual methods. Authors of |12|, addressed 
NUM problem in the context of random access in WLANs 
for stations generating either elastic or inelastic traffic. 

In this study, we focus on NUM-based rate allocation for 
streaming applications with a class of non-concave utility 
functions. Towards this, we adopt the so-called S-curve utility 
functions for streaming traffic Q, |[3l as they are shown to be 
capable of characterizing the user perceived quality for a broad 
range of multimedia streaming scenarios. In order to tackle 
the resulting non-convex NUM, we exploit transformation 
techniques to gain a strictly concave objective. However, this 



procedure yields a class of non-convex DC (difference of 
convex) constraints, referred to as reverse-convex constraints 
ifTSJ . We then deal with the non-convex transformed NUM 
using an approach called Sequential Convex Programming 
with DC constraints, abbreviated as SCP-DC, which was 
proposed in 1 13| . In this regard, SCP-DC approach tackles 
the problem with reverse-convex constraints by solving a 
series of convex problems. Then we present a distributed 
rate allocation algorithm obtained by solving the sequence of 
convex problems in an iterative manner We demonstrate that 
under mild assumptions, the proposed algorithm will converge 
to a locally optimal solution of the original NUM problem. To 
the best of our knowledge, this is the first work that addresses 
NUM with S-curve utilities with Sequential Convex Program- 
ming approach. Finally, our numerical experiments confirm 
the tractable convergence rate of our proposed algorithm and 
validate the its effectiveness in our experiment scenarios. 

The rest of this paper is organized as follows. In Section 
mi we describe the network and utility model and in Section 
Unl we establish problem formulation. Then we present our 
solution algorithm in Section HV] Numerical analysis is given 
in Section [Vl] and conclusion is drawn in Section IVIII 

II. System Model 

A. Network Model 

We consider a communication network that consists of 
a set C — {!,..., i} of unidirectional links and a set 
S = {1, . . . , 5} of sources. We denote by c = (q, Z e C) 
the link capacity vector where c/ is the capacity of link / in 
bps. We assume that each logical source s transmits at rate 
Xs ^ Xs = V^stMs], where and Ms are the minimum 
and the maximum rates, respectively. There is a fixed set of 
Unks C{s) C C that source s uses to reach its destination. 
We represent such routes using a routing matrix R G M-^^'^, 
which is defined as 



Ris 



1 if source s passes through link I 
otherwise 



Rate allocation is considered to be feasible if and only if 
the source rate vector x — {xs , s e iS) satisfies the following 
conditions 

CI. Xs ^ Xs, s G S 

C2. J2s=i RisXs <ci, I e jC. 

B. Utility Model 

In order to measure the user satisfaction degree, we use the 
well known notion of utility function. We associate an increas- 
ing and continuously differentiable function Us {xg ) with each 
source s. As mentioned in [l], multimedia applications, such 
as video streaming and VoIP, fall in the category of inelastic 
traffic and unlike elastic traffic, they are usually modeled 
by a family of non-concave utility functions referred to as 
sigmoidal-like function^ For example, previous studies 

'An increasing function f{x) is called a sigmoidal-like function, if it lias 
one inflection point xq, and f"{x) > 0, for x < xg and fix) < 0, for 
X > xq. In other words, f(x) is convex for x < xq and concave for x > xq. 
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1: S-curve Utility Corresponding to different values of Ci and 



have mostly used sigmoidal logistic function defined below as 
the utility function for inelastic traffic Q, ||9|: 



U{x) = 



1 



1 + e-a{x-[3) 



(1) 



which has the inflection point .t'"'' = /3. 

In this work, we focus on streaming applications which are 
shown to admit utility functions referred to as S-curve [31, 
111. Such utility functions can capture the perceptual video 
quality of streaming users as a function of transmission rate. 
In 1 2 1, the authors have proposed the following mathematical 
expression for this class of utility functions 

Usixs)^ , (2) 



1 



where rg is the constant rate at which video of source s 
is encoded and Xs is the average data rate received during 
transmission. Constants Cis > and C2s > 1 are some 
parameters that depend on the properties of the video sequence 
and video encoder and might be determined in an offline 
fashion for stored media streaming. 

It is easy to verify that the inflection point of the S-curve 
utility function is given by 



„infl 



C2 



1 \ "^2 = 



ClsC2 



The above equation implies that C2s > 1 results in x^"^ > 
which makes S-curve a non-concave function. Fig. [T] portrays 
some utility functions coiTesponding to different values of 
parameters Cis and C2s- 

The family of S-curve utility functions are capable of 
capturing characteristics of fine-granularity scalable, layered 
and non-scalable video streams as their special cases. For 
example, perceptual quality of FGS encoded video is smooth 
and can easily be approximated by ©. Moreover, the step- 
wise utility of SVC video streams fl4] can also be roughly 
characterized by S-curve (|2]i. Hard real-time applications such 
as traditional voice service require fixed transmission rate. For 
such services, the utility below a threshold rate would be zero. 



These applications are non-scalable and can be represented by 
a step utility function. S-curve utility function (|2|i can also 
approximate a step function as C2s ^ oo [2J. 

III. Problem Formulation 

We model the rate allocation for streaming applications 
following the framework of Network Utility Maximization 
(NUM) which was proposed as the extension to optimization 
flow control in the seminal work of Low et al. ID. The 
objective is the sum of utility functions with utilities defined 
by (|2|i and the constraints are feasibility conditions C1-C2. 
The rate allocation problem is described as follows 

Es Ris^s <ci yiec (3) 



max 



subject to 

where X denotes the Cartesian product of all rate domains 

As stated in the previous section, the S-curve utility function 
is non-concave for > 1 which makes the above 
problem non-convex. We would like to elaborate on making 
this problem convex so as to use powerful methods developed 
for convex optimization. We can make the utility function 
concave with the following change of variables: 

Substituting the above transformation in we obtain the 
transformed utility function Us {■) as 



Us (is) 



1 



e 



1 - e-Ci. 

where the transformed variable Xg belongs to 

A / ms\ I Ms 



(5) 



The transformed utility function Us is strictly concave in 
Xs, because for Cis > 0, its second derivative satisfies 

\2 



U's'iis) 



(C: 



Is 



-Cisis 



< 



1 - e-Ci= 

Rewriting the capacity constraint for link I, yields 

1 

giix) ^y^RisrsxY^ < ci; Vl e C. 



(6) 



(7) 



Unfortunately, the above capacity constraints with trans- 
formed variables do not correspond to a convex constraint 
as the (L.H.S) of © is a concave function. Indeed, the set 
Vi = {xs\gi{xs) — c; < 0} is a non-convex set, however the 
set — Vi = {xs\gi{xs) — q > 0} is a convex set. In 
optimization terminology, such a constraint is referred to as a 
reverse-convex constraint which is a special case of Difference 
of Convex (DC) constraints ITSl . lITSl . 

In order to tackle such reverse-convex constraints, we use 
the sequential convex programming algorithm with DC con- 
straints (SCP-DC) proposed in 1 13|. In this approach, the non- 
convex function that borders the range of permissible values 



for a constraint is replaced by an affine approximation to make 
the constraint convex. Using this approach, the L.H.S of each 
reverse-convex constraint gi{x) < ci is replaced by its first 
order Taylor approximation around a feasible point x, denoted 
by gi{x,x'), as follows 

gi{x,x') ^ giix') + Vgi{x'f{x - x') < ci (8) 

Since gi is differentiable, Vg exists at auxiliary variable 
x's £ Xs- It's easy to verify that gi{x,x') is affine in x 
and thereby L.H.S of (O is convex. Thus, the constraint ^ 
represents a convex constraint. For gi we get 

1 



gi{x,x') ^^Ris 



Co 



'^{Xs 



(9) 



Finally, we rewrite the NUM problem with approximated 
constraints as 



max 
subject to: 



s=l 



-Cisis 



1 



e C. 



(10) 

(11) 



gi{x,x ) < ci; 

The above problem is strictly convex (in x) since its objective 
is strictly concave because of (|6]l and its constraints are affine 
functions. 

Before proceeding to solve the above problem, it's worth 
mentioning that in case of sigmoidal logistic utility functions 
([T]i, if we define is — e°'^^^^\ we will come up with a 
convex objective with DC constraints, which can be treated by 
the aforementioned technique to obtain a convex formulation 
similar to problem (fTO]i-(fmi. Therefore, the solution procedure 
to be discussed in the next section, will be applicable to the 
case of NUM with sigmoidal logistic utility functions. 

IV. Optimal Solution 

In this section, we solve problem (fTOb- (fTTI) using dual 
methods tU, |,16J . 

A. Primal Problem 

The Lagrangian function is derived as | fT6l : 

L {x, X, ^^Us {xs) -^l^i {giix, x) - ci) (12) 

s / 

where pi is the positive Lagrange multiplier associated with 
constraint ( fTTT l for link I and fi ~ (fii,! E £). 

According to Karush-Kuhn-Tucker (KKT) theorem, the 
stationary point of the Lagrangian, i.e. the solution to 
VL (i*, x'*, /i.*) = 0, provides the unique solution to the 
problem (fTO]l-(fTni. Partial derivatives of the Lagrangian with 
respect to is is given in ( fT3T l at the top of the next page, and 
finally for the stationary point we get 



hgp" 



(14) 



where 

A* 



log 

E 



ClsC2s 



<(1- 

RisPi- 



1 — 7;— I logi' 



2s 



(15) 



ci 



Cue 



1 - e-Ci= 
1 - e-t^i= 



2s 



(13) 



It's easy to verify that the transformed source rate a;* is a 
decreasing function with respect to ^i'^ , I G C. 

We postpone finding x* to the next subsection. We will 
find /J,* by solving the dual problem associated to the primal 
problem. Towards this, we first derive the dual function, which 
is defined as the following Lagrangian maximization |fT6l : 

D{fi) = max L (^x,x , fij = L (^x* ,x'* , fij 

B. Dual Problem 

Having obtained the dual function, i.e. D (/x) = 
the dual pre 
minimization problem y_6 1 : 



L (x*,x'*,fi], the dual problem is defined as the following 



min D ifi) 



(16) 



Solving the above problem in closed form might be im- 
possible, and hence we solve it using iterative methods. As 
problem (fTO]i-(fTTTi is strictly convex, the dual function D{fj,) 
is differentiable over the open set IR+_|_ and we can benefit 
from gradient projection algorithm to solve the dual problem 
ifTTl . In this algorithm, the dual variable is iteratively updated 
in the opposite direction to \7D{fi) as follows: 

where 7 > is a sufficiently small step-size. 

Using Danskin's Theorem 1 17 1, the partial derivatives of the 
dual function are characterized as follows 

dD 



ci - gi{x, X ) 



(17) 



In an iterative setting, to find the optimal value of the aux- 
iliary variable i'*, similar to Proximal Optimization Methods 
lfT7j . we update it as follows 



;'(*+!) 



Put another way, it makes sense that, to calculate primal- 
optimal variable i*^*^^-* at iteration step t, is the best 
candidate for x along which the affine approximation (O can 
be made. 

Substituting (fTTI i into gradient projection update formula 
results in the following dual variable update 

+ 



1^1 



(18) 



where is the value of optimal transformed rate given fi 



Moreover, for rate computation at iteration t, we get 



(19) 



where 

- lo 



ClsC2s 



l-7^)logSl*) 



rs (1 - e-Ci=)^ 

Finally, by taking the inverse transformation of (|4|, for source 
rate at iterate i + 1, we get 



C2s 



(20) 



where []xs is the projection onto Xg 



V. Rate Allocation Algorithm 

A. Algorithm 

The equations obtained in the previous section for optimal 
source rate calculation, i.e. ( fT9l l and (l20l i. and dual variable 
update, i.e. dTsl ), can work together to form a distributed 
solution to problem (O. Below we have shown a concise form 
of this iterative algorithm as Algorithm 1. As we can see, 
implementation of this algorithm necessitates two mechanisms 
for information exchange between links and sources. 

1) Each link I updates its price and communicates the result 
to the corresponding sources. 

2) Each source s calculates its new rate and informs the 
links in its path. 

This type of information exchange can be carried out explic- 
itly, for example via flooding-like mechanisms as suggested in 
pSI. This is in contrast to the implicit mechanisms inherent 
in Optimization Flow Control approach H where each source 
can infer aggregate price of its path using either queueing 
delay or packet loss ratio, and each link just needs to measure 
its current flow to update its price. 

B. Convergence 

First we note that at the steady state, i.e. when 



or equivalently when 



;(*) 



F.(*+l) _ 



the approximated capacity constraints ([8]) would become 



?.(t+l) 



(21) 



and thereby convexified constraints ([Sj will become equivalent 
to DC constraints (|7]i. Therefore, if the algorithm converges to 
the steady state, the convexified constraints ([8]) will be equiv- 
alent to the original constraints of the transformed problem. 
In ifTSl . it has been proved that under mild conditions on 
the objective function, such as strict convexity, the SCP-DC 
algorithm converges to a local maximizer of the non-convex 
problem (O. Therefore, the proposed rate allocation algorithm 
will reach a local maximum of problem (O provided that 
7 is chosen sufficiently small so that the gradient projection 
algorithm will converge (TJ). 

Based on the results stated in flU, if SCP-DC algorithm 
converges, then the steady state point is a stationary point of 
the optimization problem. Put another way, the steady state 
point is a local optimal of the optimization problem. It has 
also been proved in lfT3l that for strictly convex objectives, 
the SCP-DC algorithm always converges to a KKT point of 
the original optimization problem (O. 
Algorithm 1 

Distributed Rate Control Algorithm for Streaming Traffic 
Using SCP-DC Algorithm 

Initialization 

1) Set of sources and links including the routing matrix 

2) rs,Cis,C2s for s G 5 

3) 7 and c; for I G C 

Main Loop 

Do until maxs l^s'^"'' — < e 

1) For each / G C, update link price fj,i by: 



2) For each s £ S, Xs is calculated by: 



1 



1 - 77— ] log £ i*' 

02s 



' log p- 



and p"'*' = Rig/ii^^y Then calculate Zj*^"'^'' as 



,(t+i) 



;{t + l)\ Ci 



VI. Numerical Analysis 

In this section, we investigate the performance and validity 
of the proposed rate allocation algorithm listed in the previous 
section as Algorithm 1. 

A. Scenario 1 

We first consider a simple topology with a single bottleneck 
link with capacity c = 1 Mbps. Video sequences for all 
sources are assumed to be encoded at the constant bit rate 
Ts = 256 Kbps. Sources have utility functions with parameters 
Cis = 6,Vs and (C21, . . . , C25) = (2, 4, 6, 8, 10). The 2nd 
column of Table I lists the results of the proposed rate 
allocation algorithm with step size 7 = 10^'' and the stopping 
criterion e = 0.1. 

Fig. |3(a)| and Fig. |3(b)| display the evolution of source rates 
and link price, respectively. From Fig. |3(a)[ it is observable 




Fig. 2: Network topology and flow rates 
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Fig. 3: Evolution of (a) Source rates and (b) Link price for the first 
scenario 

that the stopping criterion is met in iteration step t = \2 which 
implies the proper convergence rate of the algorithm. 

B. Scenario 2 

Now we focus on a scenario with multiple bottleneck 
links whose topology is shown in Fig. [2] with capac- 
ity vector c = (210,425,610,425,210) Kbps. Sources 
have utility functions with parameters C2S — 6,Vs and 
(Cii, . . . , C15) = (2, 4, 6, 8, 10). Similar to the previous sce- 
nario, we set rg = 256, 7 = 10^*, and e = 0.1. The 4th 
column of Table I lists the results of the proposed rate allo- 
cation algorithm. Fig. |4(a)| and Fig. |4(b)| display the evolution 
of source rates and link prices, respectively. As shown in Fig. 
|4(a)| the stopping criterion is satisfied in iteration step t = 70 
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Fig. 4: Evolution of (a) Source rates and (b) Link prices for the 
second scenario 





Scenario 1 


Scenario 2 


Source 


Algorithm 


fmincon 


Algorithm 


fmincon 


1 


118.1668 


117.9658 


209.8065 


210.0000 


2 


191.2863 


191.1745 


202.4938 


202.6043 


3 


219.4306 


219.3638 


222.3787 


227.3957 


4 


232.2970 


232.2520 


227.6564 


227.6043 


5 


239.2770 


239.2439 


197.3906 


197.3957 



TABLE L Rate allocation results 

which again demonstrates the tractable convergence rate of 
the proposed algorithm in a topology with multiple bottleneck 
links. 



C. Validation 

In order to validate the rate allocation results obtained 
above, we have also solved the problem (O by invoking 
fmincon function in Matlab 1 19|. When calling this function, 
we choose Interior-Point Method lfT6l . ifTTl as its solving 
algorithm. The results returned by fmincon for the two 
scenarios along with those obtained by our algorithm are listed 
in Table I. It is easy to confirm that the rate allocation results 
completely match those obtained from fmincon. 



VII. Conclusion 

In this paper, we addressed rate allocation for streaming 
applications with non-convex S-curve utility functions. First 
we convexified the utility functions with elementary transfor- 
mation techniques. Then, we exploited the SCP-DC approach 
ifTSl to handle the resultant reverse-convex constraints. Using 
dual method, we then proposed a distributed rate allocation 
algorithm which was shown to achieve a locally optimal 
solution of the non-convex NUM. Simulation results validated 
the tractable convergence and accuracy of the proposed rate 
allocation algorithm. As a possible direction to continue this 
research, it is promising to address rate allocation for such 
streaming applications in wireless networks. 
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